OPEN 3 ACCESS Freely available online 



•0-PLOS I o-^E 



Identification of Maize Long Non-Coding RNAs ^ 
Responsive to Drought Stress cros^rk 

Wei Zhang^ ^ Zhaoxue Han^ ^ Qingli Guo\ Yu Liu\ Yuxian Zheng\ Fangli Wu^ Weibo Jin^'^"" 

1 College of Life Sciences, Northwest A&F University, Yangling, China, 2 Institute of Bioengineering, College of Life Sciences, Zhejiang Sci-Tech University, Hangzhou, China 

Abstract 

Long non-coding RNAs (IncRNAs) represent a class of riboregulators that either directly act in long form or are processed to 
shorter miRNAs and siRNAs. Emerging evidence shows that IncRNAs participate in stress responsive regulation. In this study, 
to identify the putative maize IncRNAs responsive to drought stress, 8449 drought responsive transcripts were first 
uploaded to the Coding Potential Calculator website for classification as protein coding or non-coding RNAs, and 1 724 RNAs 
were identified as potential non-coding RNAs. A Perl script was written to screen these 1724 ncRNAs and 664 transcripts 
were ultimately identified as drought-responsive IncRNAs. Of these 664 transcripts, 126 drought-responsive IncRNAs were 
highly similar to known maize IncRNAs; the remaining 538 transcripts were considered as novel IncRNAs. Among the 664 
IncRNAs identified as drought responsive, 567 were upregulated and 97 were downregulated in drought-stressed leaves of 
maize. 8 IncRNAs were identified as miRNA precursor IncRNAs, 62 were classified as both shRNA and siRNA precursors, and 
279 were classified as siRNA precursors. The remaining 315 IncRNAs were classified as other IncRNAs that are likely to 
function as longer molecules. Among these 315 IncRNAs, 10 are identified as antisense IncRNAs and 7 could pair with 17 
CDS sequences with near-perfect matches. Finally, RT-qPCR results confirmed that all selected IncRNAs could respond to 
drought stress. These findings extend the current view on IncRNAs as ubiquitous regulators under stress conditions. 
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Introduction 

Maize {^ea mays L.) is a major cereal crop worldwide, serving as 
a major staple for both human consumption and animal feed. It 
has also become a key resource for industrial applications and 
bioenergy production. Drought is one of the major abiotic stresses 
that limit maize productivity [1]. The most effective way to 
stabilize and improve maize production under drought conditions 
is to improve the varieties in terms of drought tolerance. 
Significant differentiation of drought tolerance among maize 
genotypes implicates the hope of genetic enhancement for drought 
tolerance to improve maize [2-4] . However, breeding for drought 
tolerance is particularly challenging because of the genetic 
complexity of this trait. Drought tolerance has been well- 
documented to result from cooperative interactions among 
multiple morphological, physiological, and biochemical charac- 
ters. Different genotypes may have different responses to drought 
stress [3,5-7]. Therefore, eflicient improvement requires an in- 
depth understanding of the gene expression regulation mecha- 
nisms in response to drought stress. 

Recent genome-wide transcriptome analysis methods, such as 
tiling arrays and next generation sequencing, have revealed a large 
number of stress-responsive ncRNAs. Emerging evidence has 
revealed that ncRNAs are the major products of plant tran- 
scriptomes with significant regulatory importance [8,9]. ncRNAs 



are transcribed from intergenic regions, antisense strands of 
protein-coding genes, and pseudogenes. According to their size, 
ncRNAs are classified as small ncRNAs (sRNAs) (<40 nt) and 
long ncRNAs (IncRNAs) (>200 nt). These ncRNAs are involved 
in the transcriptional and posttranscriptional regulation of gene 
expression as well as the modulation of RNA stability and 
translation under stress conditions [10-14]. 

In contrast to sRNAs, much less is known about the large and 
diverse population of IncRNAs. This heterogeneous class of 
transcripts generally does not contain any long open reading frame 
(ORF) (no ORF >70 AA). IncRNAs can interact with proteins to 
regulate transcription, translation, or mRNA stability [15-18]. 
Furthermore, several of these long ncRNAs are precursors of 
miRNAs and siRNAs [19,20]. Similar to some miRNAs, certain 
IncRNAs are induced in various developmental processes as well 
as during abiotic stress responses in plants and animals [21-24]. In 
Caenorhabditis elegans, 25 IncRNAs are regulated in seven develop- 
mental stages and two stimulated conditions [25]. In Arabidopsis, 15 
plant IncRNAs display diverse tissue-specific expression patterns 
and/or regulation by environmental stimuli [26]. In maize, 
numerous potential IncRNAs have been identified using the full- 
length cDNA sequences of maize, and these IncRNAs may 
function to regulate expression of other genes through multiple 
RNA-mediated mechanisms [27]. However, reports on IncRNAs 
involved in drought-responsive regulation in maize are lacking. 
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In this study, we performed genome-wide screening of drought 
stress-responsive maize transcripts to isolate a collection of 
IncRNA genes expressed in maize leaves. A total of 664 putative 
transcripts were identified as drought-responsive IncRNAs. Of 
these 664 transcripts, 8 were found to be homologous to existing 
maize miRNA precursors. The remaining 656 IncRNAs were 
classified as either small RNA precursors or other ncRNAs 
through alignment with other small RNA databases. In total, 62 
IncRNAs were classified as both shRNA and siRNA precursors 
and 341 IncRNAs were classified as siRNA precursors, including 
62 shRNAs. The remaining 315 IncRNAs were classified as other 
IncRNAs that are likely to function as longer molecules. Among 
these 315 IncRNAs, 1 0 were identified as antisense IncRNAs and 7 
could pair with 17 CDS with near-perfect matches. qRT-PCR 
results confirmed that all selected IncRNAs could respond to 
drought stress. These findings extend the current view on 
IncRNAs as ubiquitous regulators under stress conditions. 

Materials and Methods 

Data sets 

Maize genomic sequences (B73_RefGen_v2), genome annota- 
tion-filtered site file (ZmB73_5b_FGS), and the sequences of 
63540 CDSs were downloaded from the ftp site: ftp: //ftp. 
maizesequence.org/pub/maize/release-5b/. Small RNA sequenc- 
es were downloaded as FASTA files from NCBI under the series 
GSE15286 (http:/ /www.ncbi.nlm.nih.gov/geo/query/acc. 

cgi?acc = GSE 15286) [28]. Files for both roots and shoots for 
each group (known_miRNA, shRNA, and siRNA) were combined 
into a single FASTA file. All of the miRNA precursor sequences 
were downloaded from miRBase (Version 20) [29] at the website: 
http://www.mrrbase.org/. BLAST+ executables, version 2.2.23, 
was downloaded from the NCBI website for local use. 

Bioinformatics analysis 

AH sequenced reads from each sample were aligned to the maize 
reference genome (B73_RefGen_v2) using Tophat [30]. The 
transcriptome of each sample was assembled separately using 
Cuflfiinks [31]. Then the assembled transcript isoforms were 
compared with the maize genome annotation-filtered site (FGS), 
and producted a non-redundant genome features annotation file 
(Generic Feature Format, gff) [32]. Only reads uniquely mapped 
to the non-redundant gff annotated site were kept for expression 
analysis. The expression level of every isoform was calculated by 
using RPKM method (Reads per kHobase transcriptome per 
million mapped reads) [33]. The sequences of drought stress 
responsive transcripts, which were 2-fold up- or downregulated 
with False Discovery Rates (FDRs) [34] less than 0.001, were 
retrieved using a Perl script, and then were uploaded to the 
Coding Potential Calculator (CPC) website for classification as 
protein-coding or ncRNA candidates. Finally, a Perl script was 
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Figure 1. The flowchart for predicting non-coding RNA 
responsive to drought stress. 

doi:10.1371/journal.pone.0098958.g001 

written to screen these ncRNAs for RNA length >200 nt and 
ORF <80 AA, and drought-responsive IncRNA sequences were 
extracted. The identification flowchart of IncRNAs is shown in 
Figure 1. 

To identify the drought-responsive IncRNAs as precursors of 
known miRNAs, the miRNA precursors were aligned to IncRNAs 
using BLAST. The IncRNAs matched the miRNA precursors with 
>99% identification accuracy; those IncRNAs with >90% 
coverage were finally classified as miRNA precursors. 

To identify drought-responsive IncRNAs as precursors of 
shRNAs or siRNAs, the shRNA and siRNA libraries were 
respectively mapped to the IncRNAs using Bowtie without 
mismatches [35]. IncRNAs homologous to the shRNA reads were 
classified as shRNA precursors while those homologous to the 
siRNA reads were classified as siRNA precursors. 

Plant growth and drought stress treatment 

Maize [^ea mays cv B73) seeds were soaked in deionized water 
for 1 2 h and then placed on a sheet of moist filter paper in a Petri 



Table 1. Identification and expression analysis of IncRNAs responding to drought stress. 






Up-regulatlon 


Down-regulation 


total 


Intergenlc IncRNAS 


344 


77 


421 


Intragenic IncRNAs 


189 


16 


205 


Antisense IncRNAs 


21 


2 


23 


Overlapping IncRNAs 


13 


2 


15 


Total 


567 


97 


664 
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Figure 2. Quantitative analysis of the precursor levels of 5 miRNAs in the leaves using RT-qPCR at 0, 5 and 1 0 h. GAPDH RNA was used 
as the internal control. Error bars Indicate SD obtained from three biological repeats. 
doi:1 0.1 371/journal.pone.0098958.g002 



dish. These seeds were germinated at 28°C for 3 d. Germinated 
seeds were transferred to a floating foam sheet in the hydroponic 
boxes (40 x20 x 12cm'^) containing continuously aerated water in a 
growth chamber (28°C day/26°C night, 16 h photoperiod, 30%- 
50% relative humidity) for about 1 wk. The seedlings were then 
cultivated in 1/2 Hoagland solution [2.5 mM KNO3, 2.5 mM 
Ca(N03)2-4 H2O, 1 mM MgS04-7 H2O, 0.5 mM KH2PO4, 
50 ^iM Fe-EDTA, 7.5 ^iM H3BO3, 2.5 ^iM (NH4)6Mo7024, 
1.25 |J,M MnCl2, 1 |iM ZnS04, and 0.5 |iM CUSO4, pH 6.0] 
for about another week. To induce expression of target genes, 
seedlings at the three-leaf stage were subjected to drought stress 
treatment. Drought treatment was carried out by submerging the 
roots of the plants in 1/2 Hoagland solution with 16% (w/v) 
polyethylene glycol (MW 8000) for different periods. The shoot 
and root tissues of the control and stressed seedlings were 
harvested at 3 time points (0 h, 5 h, and 10 h) after treatment in 
three biological replicates. These samples were immediately frozen 
in liquid nitrogen and stored at — 70°C for studies on dynamic 
expression changes in the transcripts. 



RNA isolation and quantitative real-time PCR analysis 

Total RNAs were extracted from leaf tissues using a TRIzol 
reagent (Invitrogen, Carlsbad, CA, USA), followed by RNase-free 
DNase treatment (Takara, Dalian, China). The RNA concentra- 
tions were quantified by a NanoDrop ND-1000 spectrophotom- 
eter. 

The expression profiles of drought-responsive IncRNAs were 
assayed by Reverse Transcription (RT) quantitative PCR (qPCR). 
500 ng of total RNAs were used for initiating the RT and the RT 
product was used as template for qPCR using the IncRNAs 
specific primers. The RT reactions were performed using M-MLV 
Reverse Transcriptase (Takara, Japan) according to the supplier's 
protocol. Primers were then added to perform PCR. The GAPDH 
was used as the internal control for RT-qPCR. All the oligos used 
in this study were listed in Table S 1 . 

SYBR Green PCR was performed as per the manufacturer's 
instructions (Takara, Japan). Briefly, 2 |J,1 of cDNA template was 
added to 12.5 |xl of 2x SYBR Green PCR master mix (Takara), 
1 |J,M concentration of each primer, and ddH20 to a final volume of 
25 |xl. The reactions were ampHGed for 10 s at 95°C, followed by 40 
cycles of 95°C for 10 s and 60°C for 30 s. AU reactions were 



Table 2. The IncRNAs were identified as miRNA precursors. 



Gene_ID 


maize mlRNA precursor 


% coverage 


% Identity 


Up or Down On RAN-Seq 


TCONS_00012662 


MIR167j 


100 


100 


Up 


TCONS_00046477 


MIR169d 


100 


100 


Down 


TCONS_00042984 


MIR169h 


100 


100 


Up 


GRMZM2G420571_T01 


MIR172C 


100 


100 


Up 


TCONS_00044116 


M[R399b 


100 


100 


Down 


TCONS_00024310 


MIR399e 


100 


100 


Down 


TCONS_00024309 


MIR399e 


100 


100 


Down 


TCONS_00054544 


MIR827 


100 


100 


Down 
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TCONS_00012703 

t0613221_l ftTGCTCCffPGATCTTTTTOGGT , 

tl592118_l ACGTATAGGTAATCATCTCGAAA 

t0747e08_l CGTATAOGTJiATCATCTCaAAAGA 

t0067360_2 QATQTAGACCTCATAGTGATAGGC 

t0072913_3 ATGTAGACCTCATAGTGATAGOCA 

t0204396_l TGTAGACCTCATAGTGATAGGCAT 

t0344ia9 1 TOTAQACCTCATAaTOftTAGGCAT 

t02B6653_l TTTCAAATGCGGCCTAAAGCACGG 

t095e937_l TTCAAATQCGGCCTAAAGCACGGT 

tl215302 1 AAATGCGGCCTAAAGCAC6GTCTC 

tl414078_l ATOOCJGCCrAAAaCROQOTCTCTQ 

t0249345_l ATAAAOAATTTAacaCAaCCAGT. . 



Figure 3. The TCONS_00012703 were mapped with multiple siRNA reads. 

doi:1 0.1 371/journal.pone.0098958.g003 



performed in triplicate, and the controls (no template and no RT) 
were included for each gene. The relative expression levels were 
calculated using the formula: Ratio = (iWget) ^'-^P^^g^j^'^"'™'"'"**™''"''/ 

(i^ef) '^Gp_^/'^ontrol-treatment)|-3g-|_ 



Results 

Identification and expression characterization of drought 
stress-responsive long ncRNAs in maize 

A total of 8449 drought stress responsive transcripts that were 2- 
fold up- or downregulated with FDRs less than 0.001 were 
obtained from RNA-seq data (Unpublished data from Zhaoxue 
Han's lab). In thi.s study, to identify putative maize IncRNAs 
responsive to drought stress, the 8449 RNA sequences were first 
uploaded to the CPC website for classification as protein coding or 
non-coding RNAs; 1724 RNAs were identified as potential non- 
coding RNAs (Data SI). A Perl script was then written to screen 
these 1724 ncRNAs with the RNA length >2()0 nt and ORF <8() 
AA; a total of 664 transcripts were finally identified as drought- 
responsive IncRNAs including 205 intragenic IncRNAs, 421 
intergenic IncRNAs, 23 antisense and 15 IncRNAs overlapping 
with parts of inter- and intragenic sequences (Figure 1, Data S2, 
GenBank accession number: KJ682450-KJ682628, KJ731849- 
KJ732333). The expression profiles of these 664 IncRNAs 
obtained between control and drought stress samples were 
extracted from RNA-seq data (Table S2). Among the 664 
IncRNAs, 567 IncRNAs containing 344 of intergenic, 189 of 
intragenic, 21 of antisense and 13 of overlapping IncRNAs were 
upregulated in drought-stress leaves. The remaining 97 were 
downregulated including 77 of intergenic, 16 of intragenic, 2 
antisense and 2 overlapping IncRNAs (Table 1, Table S2). 
Moreover, 1 26 of the 664 drought-responsive IncRNAs showed 
high homology (>94% identity and >80% coverage) to known 
IncRNAs identified by Boerner et al. [27] and were thus 
considered as known IncRNAs (Table S3). 

ncRNA transcripts correspond to miRNA precursors 

By aligning miRNA precursors to the 664 long ncRNAs, we 
identified 8 IncRNAs as 7 known maize miRNA precursors of 
miR167j, miR169d, miR169h, miR172c, miR399b, miR399e and 
miR827 (Table 2). Among the 8 IncRNAs, 
GRMZM2G420571_T01 was an intragenic IncRNA and the 
TCONS_00054544 was an overlapping IncRNAs, and the 
remaining 6 were intergenic IncRNAs. To prove that these 
miRNA precursor IncRNAs could respond to drought stress, 
quantitati\'e RT-PCT (qRT-PCR) was employed to detect their 
expression in both drought-stressed and control samples. The 
precursor sequences of miR169d and miR169h are high similarity 
and the expressions of them were opposite on RNA-seq data. 
Therefore, the expression detections for the other 5 miRNA 
precursor IncRNAs were performed and the results show the 
obviously upregulated in stressed leaves compared with control 
leaves (Figure 2). Moreover, only 2 IncRNAs, namely the 



precursors of miR167j and miR172c, exhibited consistent 
expression changes compared with the RNA-seq data. The 
remaining 3 miRNA precursors exhibited upregulation in 
drought-stressed leaves as determined by qRT-PCR but showed 
downregulation in comparison with RNA-seq (Table 2). There- 
fore, these results suggest that miR167j and miR172c precursor 
IncRNAs could be responsive to drought stress. 

Characterization of ncRNAs for shRNA and siRNA 

To identily drought-responsive ncRNAs as precursors of 
shRNAs or siRNAs, small RNA libraries of shRNAs and siRNAs 
were respectively mapped to the remaining 656 IncRNAs using 
Bowtie without mismatches. The results show that 62 IncRNAs are 
homologous to numerous shRNA reads and may hence be 
classified as shRNA precursors (Table S4). About 341 IncRNAs 
were homologous to numerous siRNA reads and hence classified 
as siRNA precursors (Table S5). In addition, we found that all of 
62 shRNA precursors are homologous to siRNA library reads and 
may hence be classified as siRNA precursors. The ncRNAs are 
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Figure 4. Quantitative analysis of 6 siRNA precursor IncRNAs in 
the leaves and roots using RT-qPCR at 0, 5 and 10 h. The gray 
bar denotes leaf samples and the white bar denotes root samples. 
GAPDH RNA was used as the internal control. Error bars indicate SD 
obtained from three biological repeats. 
doi:l 0.1 371/journal.pone.0098958.g004 
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Table 3. Classification of 664 drought-responsive IncRNAs. 







intergenic 


Intragenic 


antisense 


overlapping 


Total 


miRNA 


6 


1 


0 


1 


8 


shRNA/siRNA 


52 


5 


4 


1 


62 


siRNA 


196 


67 


9 


7 


279 


Other IncRNAs 


167 


132 


10 


6 


315 


Total 


421 


205 


23 


15 


664 
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homologous to multiple types of small RNAs, which suggests 
overlap between sRNA biogenesis pathways at numerous loci. In 
plants, evidence of overlapping epigenetic regulatory pathways 
exists [37,38]. 

Based on the established mechanisms of biogenesis, siRNAs are 
expected to derive from a longer molecule cleaved into multiple 
sRNAs through specific endonuclease activity. To determine if 
multiple sRNAs are derived from siRNA precursor ncRNAs, 341 
siRNA precursor ncRNAs were mapped to siRNA reads and 
subjected to further analysis. Results show that most of the siRNA 
precursor IncRNAs produced more than one sRNA (Table S5). 
For instance. Figure 3 shows that TCONS_00012703 is homol- 
ogous to multiple siRNA reads. In addition, the RT-qPCR result 
of TCONS_00012703 RNA demonstrates that the expression 
level of the IncRNAs is relatively low in both non-stressed leaves 
and roots but accumulates to high levels after drought stress in 
both plant tissues (Figure 4A). 

To validate these drought-responsive IncRNAs, qRT-PCR was 
performed to detect the expression profiles of 5 siRNA precursor 
IncRNAs randomly selected from the 279 IncRNAs at 3 time 
points of 0 h, 5 h and 1 0 h after drought treatment. Similar to the 
RNA-seq data (Table S2), the RT-qPCR results show that 4 
siRNA precursor IncRNAs are upregialated in stressed leaves 
(Figure 4 B-E) except for a root-specific IncRNA of 
TCONS_00056395 (Figure 4F), and exhibit maximum expression 
in leaves stressed for 5 h. In addition, TCONS_00082 1 74, which 
was only expressed in leaves, was not detected in roots, which 
suggests leaf-specific expression (Figure 4E). 
GRMZM2G088590_T04 which was downregulated in root 
stressed for 5 h, was upregulated in roots stressed for 10 h 
(Figure 4D). The remaining 2 IncRNAs, TCONS_00037470 and 
TCONS_00012768, were downregulated in drought-stressed 
roots (Figure 4B and C). These findings confirm that these siRNA 
precursor IncRNAs are responsive to drought stress in maize 
leaves and roots. Moreover, based on the differential expression 
profiles observed between leaves and roots, these IncRNAs may 
participate in gene expression regulation in different pathways 
within leaves and roots. 

Characterization of the remaining 315 IncRNAs 

Among the 664 drought-responsive IncRNAs, the 3 1 5 IncRNAs 
that were not classified as small RNA precursor IncRNAs, 
including miRNAs, shRNAs and siRNAs precursors, were finally 
classified as other IncRNAs that are likely to function as longer 
molecules (Data S3). Among them, 167 were intergenic IncRNAs, 
132 were intragenic IncRNAs, 10 were antisence transcripts and 
the remaining 6 were overlapped with parts of intragenic and 
intergenic regions (Table 3). To determine whether or not these 
IncRNAs could regulate the expression of protein coding genes as 



long molecules, excepted for 1 0 antisense IncRNAs, the remaining 
305 IncRNAs were tested for homology to CDS sequences using 
BLAST with Plus/Minus patterns with >90% identity and >90% 
coverage. The results demonstrate that 7 IncRNAs could pair with 
17 CDS sequences with near-perfect matches (Table 4), which 
suggests that these 7 IncRNAs might regulate the expression of 1 7 
proteins by inducing transcriptional or post-transcriptional gene 
silencing. 

IncRNAs have been found to exhibit a wide range of functions. 
The transcription of certain IncRNAs is highly tissue- and 
temporal-specific, and their expressions may be responsive to 
certain stimuli. To analyze the tissue-specific expression of these 
315 IncRNAs, the expression data of 70061 maize transcripts in 9 
tissues (10DAP_Whole_Seed, Rl_Anthers, Rl_Innermost_Husk, 
Rl_Pre-poUination_Cob, Rl_Silks, Vl_GH_Primary_Root, 
Vl_Pooled_Leaves, V3_Topmost_Leaf and V4_Stem_and_SAM) 
were downloaded from the maizeGDB website, http://maizegdb. 
org/, and reversed to the raw expression values. The raw 
expression values of the 9 tissues were extracted based on GenelD. 
Ultimately, only 57 of the 315 drought-responsive IncRNAs could 
be obtained from the tissue expression data. According to the 
expression profiles of these 57 IncRNAs in 9 tissues, most IncRNAs 
exhibited tissue-specific expression (Figure 5, Table S6). 

Four IncRNAs, GRMZM2G574383_T01, TCONS_000 12690, 
TCONS_00007700 and TCONS_00000649 were randomly 
selected from the 315 IncRNAs for expression analysis in maize 
leaves and roots at 3 time points of 0, 5, and 10 h after drought 
treatment using qRT-PCR. Similar to the RNA-seq data (Table 
S2), RT-qPCR results showed that the 4 IncRNAs were obviously 
upregxilated in stressed leaves (Figure 6). TCONS_00000649 
IncRNA exhibited leaf-specific expression, whereas the expression 
of the GRMZM2G574383_T01 was relatively low in leaves but 
very high in roots, suggesting a root-specific IncRNA. Further- 
more, TCONS_00007700, which was upregulated in stressed 
leaves, was downregulated in stressed roots (Figure 6). 

Discussions 

LncRNAs are a type of novel molecule with important functions 
in a wide range of biological processes, including developmental 
regulations and stress responses; nevertheless the detailed mech- 
anisms involved in this biological processes remain largely 
unknown [39]. In this study, to identify the putative maize 
IncRNAs responsive to drought stress, a total of 664 out of 8449 
differentially expressed transcripts in maize were identified as 
drought-responsive IncRNAs, including 126 known and 538 novel 
IncRNAs. In plants, small RNAs are an important class of 
noncoding RNAs for the regulation of gene expression at the 
transcriptional or posttranscriptional level. These sRNAs can 
originate from longer precursors that are processed by endonu- 
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Figure 5. The tissue specific expression of 57 IncRNAs. 

doi:1 0.1 371 /journal.pone.0098958.g005 

cleases, such as those identified by Dicer [40] . In this study, 8 out 
of the 664 drought-responsive IncRNAs were identified as 
precursors of miR167j, miR169d, miR169h, miR172c, miR399b, 
miR399e and miR827. RT-qPCR also confirmed that miR167j 
and miRI72c precursor IncRNAs participate in maize drouglit- 
stress responses. These miRNAs have been reported as drouglit- 
responsive miRNAs in rice or maize in previous studies [41,42]. 

Among the 664 drouglit-responsive IncRNAs determined, 315 
IncRNAs including 10 antisenses and 305 non-antisenses were 
ultimately identified as other IncRNAs that are likely to function as 
longer molecules. To understand their function, the 305 non- 
antisense IncRNAs were tested for homology with CDS sequences. 
Results show that 7 IncRNAs could pair with 1 7 CDS sequences, 
which suggests that these 7 IncRNAs regulate the expressions of 1 7 
proteins by inducing transcriptional or post-transcriptional gene 
silencing. Antisense transcripts can mediate gene silencing via 
transcriptional or post-transcriptional mechanisms, with the latter 
involving mRNA degradation. Antisense-mediated mRNA degra- 
dation likely generates siRNAs, and these siRNAs were not 
detected in the available datasets. The IncRNA candidates may 
produce siRNAs, but these siRNAs are not represented in the 
datasets used to create the databases for this study. 

Numerous IncRNAs have been reported to participate in 
responses to a wide variety of stresses, including biotic or abiotic 
stresses. In Arabidopsis, the expression of 1832 IncRNAs has been 
changed after 2 h and/ or 1 0 h of drought, cold, high-salt, and/ or 
abscisic acid (ABA) treatments [43]. In wheat, 125 IncRNAs were 
also identified as powdery mildew-responsive IncRNAs and heat 
stress-responsive IncRNAs [44] . In this study, some IncRNAs were 
selected through their expression in both drought-stressed and 
control samples using RT-qPCT. Results demonstrate that all of 



the selected IncRNAs show apparent changes in their expression 
profiles in leaves stressed for 5 h or 10 h. In addition, several 
IncRNAs exhibited leaf-specific or root-specific expression. These 
results prove a previous finding that IncRNAs respond to abiotic or 
biotic stress via tissue-specific expression [25,26]. 
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Figure 6. Quantitative analysis of 4 other IncRNAs in the leaves 
and roots using RT-qPCR at 0, 5 and 10 h. The gray bar denotes 
tfie leaf samples and tfie wfiite bar denotes root samples. GAPDH RNA 
was used as tfie internal control. Error bars indicate SD obtained from 
three biological repeats. 
doi:10.1371/journal.pone.0098958.g006 
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